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Abstract. Electronic properties, band width, band gap and van Hove singu- 
larities, of (3,0), (4,0) and (9,0) zigzag nanotubes are comparatively investi- 
gated in the Harigaya's model and a toy model including the contributions of 
bonds of different types to the SSH Hamiltonian differently. Optical transition 
frequencies are calculated. In this way an experimental correlation between 
the two models is achieved. 



1. Introduction 

The last two decades have seen an explosive growth of the investigations on 
nanoscale materials. Carbon nanotubes have taken a rather large part of this effort 
by blazing a trail in experimental and theoretical studies. This is due to their unique 
chemical, electrical and mechanical properties because of the one-dimensional (ID) 
confinement of their electronic states, resulting in the so-called van Hove singu- 
larities (vHSs) in the density of states (DOS). With their engrossing properties 
they have enormous applications in very different areas, ranging from electron- 
ics to biotechnology. Research on controlled synthesis of nanotubes to progress 
the standards of appHcations is still a grand challenge. Recent progresses may 
make single-wall carbon nanotubes (SWCNTs) to play an important role in future 
nanometer scale integrated photonic devices, in quantum optics and in biological 
sensing [1]. 

SWCNTs can be thought of as graphite cylinders formed by connecting two crys- 
tallographically equivalent sites in a graphite sheet. A graphite sheet is a 7r-system 
based on sp^ hybridized carbon atoms consisting of strongly bound tr-electrons and 
a zero-gap semiconductor in the sense that the conduction and valance bands con- 
sisting TT states cross at K and K' points of Brilloin zone. As derived from the 
graphite sheet by the zone folding method, a (n, m) type tube is metallic, possess- 
ing nonzero DOS in the region near Fermi energy, if n — to is a multiple of 3 since 
electronic band crosses the Fermi level at two-thirds of the way to the zone edge. 
These are correct as far as we remain true to zone folding approximation [2]. 

The tight - binding approximation has played an important role as one of the 
most commonly used calculation method for quasi - one - dimensional and higher- 
dimensional systems. The Su-Schrieffer-Heeger (SSH) model acquiesces this ap- 
proximation and is commonly used to study the stability and electronic structure 
of nanotubes. In this model, hops of the tt electrons between neighboring car- 
bon atoms are studied and the local electron-phonon interaction is treated in an 
adiabatic approximation [3,4]. 
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The SSH model treats the roles of all bonds, tilt or right, on an equal foot- 
ing. However, chainlike structure of graphene and hence of nanotubes brings in 
mind that the role of tilt and right bonds in the lattice and electronic structure 
of nanotubes might be different. As a matter of fact Cabria et al. [5] during an 
investigation of the stability of narrow zigzag carbon nanotubes in the context of 
local density functional calculations found that the axial bonds becoming shorter 
and more sp^-like in character, while the other bonds are becoming longer and more 
sp^-like in character. 

Inspired by these, we found it worthwhile to modify the SSH Hamiltonian used 
by Harigaya for studying the stability and energy gap properties of nanotubes [6-10] 
by considering the contributions of bonds of different types to the SSH Hamiltonian 
differently at the very beginning. We introduced this modified version of Harigaya's 
model and its appHcation to armchair type nanotubes in [11] and some other pre- 
liminary results elsewhere [12]. Concerning the (3,0) zigzag nanotube, we observed 
that the tiny energy gap appearing in Harigaya's model was lost, a result consistent 
with the metaUicity of this tube [13]. In this work, we continue to investigate the 
electronic structure of zigzag type nanotubes. 



2. Model and Calculation 

The aforementioned modified version of Harigaya's model is based on the SSH 
Hamiltonian which includes different roles of the tilt and right bonds. This is 
achieved by taking the hopping integral of the ideal undimerized system, the electron- 
phonon coupling constants and the spring constants different for the tilt and right 
bonds. In this way the introduction of two Lagrange multipliers and hence two con- 
straints become unavoidable. The separation of the constraint into two constraints, 
vanishing of the sum of right bond distortions and vanishing of the sum of tilt bond 
distortions, provides more freedom for lattice relaxations. 

The total energy of the system is read as 



i,(T i i 

where €j,aS are the eigenvalues of the model Hamiltonian. The self- consistent 
equation for the lattice is extracted as 



(2) 



= 27*^ 



Above Bs arc the eigenvectors of the Hamiltonian and 7*''' = [q;*^'']^/k*'''. The 
prime indicates the sum over the occupied states of electrons and iVb is the total 
number of tt- bonds. It is 1.5 times of the number of sites N. Et and vf^'s are 
numerically calculated in the iteration method by the adiabatic approximation for 
phonons and the calculation was kept further to have converged once vl'^ varied 
less than 10~^. 
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3. Results and Discussions 

The determination of energy eigenvalues, ei^^s and energy eigenstates, Bs by 
solving Schrodinger's equations necessitates the knowledge of the matrix represen- 
tation of model's Hamiltonian. It is quite clear that this matrix representation 
depends strictly on the geometry, i.e., whether one plans to do the calculations for 
a carbon sheet, for an open ended nanotube or a nanotube with periodic boundaries, 
since the Hamiltonian is different for different geometries. 

To be able to calculate the electronic band structure properties of different ge- 
ometries all together we preferred to write the matrix representation of electronic 
part of our Hamiltonian in a compact manner as follows 



(3) 
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Here, A, B, C and D stand for the hopping integrals between the nearest - 
neighbor interaction sites in the same row (tilt bonds) ; between the nearest - neighbor 
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interaction sites in different rows (right bonds), which causes the sheet structure; 
between the first and last sites in the same row, which creates zigzag nanotube, and 
between the first and the last rows, which create the periodic boundaries, respec- 
tively. All the remaining terms indicated by dots are zero. /3s, Pn and /3t represent 
the evolution parameters for carbon sheet, open ended nanotube and nanotube 
with periodic boundaries. Each beta parameter varies from to 1. In the case of 
armchair nanotubes the Bs and Ds exchange. 

Because of our toy model gives a better result concerning the electronic band 
structure of (3,0) nanotube, compared to the Harigaya's model [13], we want to 
check its results from the point of view of the other electronic structure properties 
of this nanotube and also from the point of view of electronic structure properties 
of a large radius conducting (9,0) nanotube together with a semiconducting (4,0) 
nanotube. We perform the calculations with N = 2n and K up to 100 for all tubes. 
For (3,0), the total number of C-sites increases up to 600 for K = 100 and the 
numbers of tilt and right bonds reach up to 600 and 297. These numbers are 800, 
800 and 396 for (4,0) while they are 1800, 1800 and 891 for (9,0). In the periodic 
boundaries case, the numbers of right bonds keep themselves while the numbers of 
tilt bonds increase 3, 4 and 9 more, respectively. Of course, it should be mentioned 
that for zigzag tubes the integer n is related to the length while the integer K is 
related to the diameter of tube. We chose the numerical set of t^^ , a*''' and 
parameters as 2.5 eV, 6.31 eV/A and 49.7 eV/A^, respectively [10]. 

The energy eigenvalues when (3,0), (4,0) and (9,0) nanostructures (graphene, 
open ended nanotubes and nanotubes possessing periodic boundaries) evolve from 
trans-PA chains, both in the framework of our toy model and in that of Harigaya's 
model, are obtained numerically by taking to and a values the same for tilt and 
right bonds. From these eigenvalues we extracted the band width and band gap 
variations together with the one-dimensional vHSs which correspond to logarithmic 
discontinuities in DOS. 

First, we report band width, i.e.; the difference between the highest and lowest 
energy eigenvalues. Figs. 1(a), (b) and (c) show band width for three systems 
comparatively in both models. For (3,0) the band width naturally starts with 8.73 
eV, the band width value of trans-PA chain having the suitable length. According 
to the Harigaya's model, at the early stages of carbon sheet formation band width 
rises steeply and it then enters a sequence of rising and staying constant stages with 
small periods. But, according to our toy model it regularly increases during the 
whole sheet formation process. However, at the end of sheet formation process, band 
width reaches the same value, 13.80 eV for both models. During the completion 
of open ended nanotube both models give the same behavior for band width. It 
reaches 15 eV. According to both models, the band width stays almost constant 
at the formation stage of nanotube with periodic boundary conditions. If only the 
final value of band width would be considered, there will not be any difference 
between the two models. For (4,0) the band width starts with 9.16 eV, again the 
band width of trans-PA of length yielding the (4,0) structures, and shows almost 
the same behavior as the band width of (3,0). For (9,0) the same appears with 9.75 
eV, the band with of trans-PA of this length. 

Next, we show band gap, i.e., the difference between HOMO and LUMO. For 
(3,0) the gap initiate with 3.39 eV, the gap value of trans-PA. As graphene evolves, 
the gap goes down almost in a parallel way for both models to 0. This is consistent of 
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Figure 1. Band with evolution of (3,0), (4,0) and (9,0) structures 
both in the Harigaya's model and in the toy model. As /3s varies 
both /?N and /3t are kept fixed at zero; as /3n varies /3s is kept fixed 
at 1 and /3t is kept fixed at zero and as /3t varies both /3s and 
/3t are kept fixed at 1. Solid and dotted lines correspond to the 
variations in the Harigaya's model and in the toy model, respec- 
tively. 



graphene being a zero-gap semiconductor whose electronic structure near the Fermi 
energy is given by an occupied tt band and an empty tt* band. It is quite clear from 
Fig. 2(a) that in the Harigaya's model the gap drops to zero for = 1-63 eV while 
it drops to zero for ^q"^ = 1.38 eV in the toy model. This is due to the looseness 
of the constraint conditions compared to those in Harigaya's model. At the latest 
stages of graphene formation and during the formation of open ended nanotube, 
the gap vanishes and with the beginning of nanotube with periodic boundaries it 
jumps up to 0.39 eV and falls to 0.33 eV at the completion when Harigaya's model is 
used. Hence, according to this model we have a nonzero gap for a (3,0) nanotube. 
SSH model does not include bond angle and bond length variations giving rise 
to finite curvature effects in nanotubes, detected in several experiments as small 
energy gap near the Fermi level [14]. Therefore, the tiny gap near the Fermi level 
in (3,0) appearing according to Harigaya's model can not be interpreted as due to 
curvature effects. For this an extended form of tight binding calculations to include 
curvature effects had to be used [15]. But according to our toy model, the gap starts 
and ends with 0.08 eV and 0.07 eV during the formation stage of (3,0) nanotube 
with periodic boundary conditions (Fig. 2(a)). For (4,0) the initial gap value is 2.89 
eV, it follows almost the same decrease up to zero and then with the beginning of 
nanotube with periodic boundaries jumps to relatively close values, 2.23 eV in the 
Harigaya's model and 2.07 in the toy model (Fig. 2(b)). Peng et al. [16] assume 
that this SWCNT is a semiconductor, which was predicted to have an energy gap 
of about 2.5 eV based on the graphene sheet model. On the other hand, some 
authors [5] claim that (4,0) is metallic and this discrepancy arises from the neglect 
of curvature effects in the graphene sheet model. As for (9,0), the gap starts from 
an even smaller value, 1.91 eV and then falls to zero at to = 0.88 eV and at — 75 
eV, respectively. At the completion of nanotube with periodic boundaries the gap 
goes to zero (Fig. 2(c)). 

We will report on calculations of one-dimensional van Hove singularities, i.e., 
logarithmic singularities: These singularities may be important in experimental 
correlations of models of nanotubes because optical absorption or emission rates 
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Figure 2. Band gap evolution of (3,0), (4,0) and (9,0) structures 
both in the Harigaya's model and in the toy model. As /3s varies 
both /3n and /?t are kept fixed at zero; as /3n varies /3s is kept fixed 
at 1 and /3t is kept fixed at zero and as /?t varies both /3s and 
/3t are kept fixed at 1. Solid and dotted lines correspond to the 
variations in the Harigaya's model and in the toy model, respec- 
tively. 



in nanotubes is related primarily to the electronic states at the vHSs [17]. In this 
respect, the electronic transition energies En, i — 1,2,3,..., that is the energy 
separations between i"* vHSs in the DOS of the conduction and valence bands play 
the fundamental role since the measured absorption peaks occur at these transition 
energies [18,19]. 

The recent progress in measuring the absorption peaks of single SWCNT instead 
of SWCNTs bundle is encouraging [1] . In this way effects arising from ensemble av- 
eraging and effects introduced by a composite medium are eliminated. Besides this, 
combined fiuorescence and Raman scattering measurements on the single molecule 
level open up new ways to characterize SWCNTs [1]. These are among nondestruc- 
tive characterization techniques of SWCNTs, the only handicap is the effects on 
electronic properties due to inhomogeneities in their local environments. 

Sharp vHSs occur in carbon nanotubes with diameters less than 2 nm [6]. All of 
our three SWCNTs have this property, 0.235, 0.313 nm and 0.704 nm, respectively. 
Figs. 3, 4 and 5 respectively depict the DOS of graphene, open ended SWCNTs and 
nanotubes possessing periodic boundaries according to the Harigaya's model and 
the toy model for (3,0), (4,0) and (9,0) structures. We calculated the energies at 
which these singularities appear. In the toy model the number of vHSs is higher. 
Gap is the interval between the first vHSs in the valance and conduction bands. For 
(3,0), the gap value obtained in this way is consistent with the gap value obtained 
from Fig. 2(b) both in the Harigaya's model and the toy model. (4,0) has an 
energy gap of about 2.5 eV based on the graphene sheet model [16]. Our result is 
consistent with this. For (9,0), we obtain an extra van Hove singularity very close 
to Fermi energy according to both models. This is reasonable since (9,0) nanotubes 
are conducting. The next three singularities are coinciding with those of Hartschuh 
et al. [1]. 

Transition frequencies uuu = Eu/h cm~^ in the optical spectral range for (3,0), 
(4,0) and (9,0) open ended nanotubes and nanotubes with periodic boundaries are 
selected from the calculated transition energies and are shown in Table 1. 
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Figure 3. DOS evolution of (3,0 ) structures for K = 60. (a) for 
the Harigaya's model and (b) for the toy model. 
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Figure 4. Same as Fig. 3 but for (4,0 ). 



4. Conclusion 

Electronic properties, band with; band gap and vHSs of three zigzag nanotubes, 
(3,0) (small radius-conducting), (4,0) (semiconducting) and (9,0) (large radius- 
conducting) are comparatively studied in a toy model and in the Harigaya's model. 
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Figure 5. Same as Fig. 3 and Fig. 4 but for (9,0 ). 



The toy model is an extension of the Harigaya's model, it includes the contributions 
of bonds of different types to the SSH Hamiltonian differently. Both models give 
the same band width. In the (3,0) case, the band gap appearing according to the 
Harigaya's model disappears when the toy model is used for periodic conditions. 
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Tubes 


Models 




Frequencies (cm ^) 




Harigaya 


Pn = 


1 


25249.3 












(3,0) 


ioy 


PT = 

Pn = 


1 

1 


ZZLboA 
19988.1 


2o887.7 
26858.3 












M 


Pt = 


1 


20910.7 


23083.6 


26941.8 










Harigaya 




1 


19881.8 


24102.7 


28094 








(4,0) 


II 

Toy 

M 


/3t = 
Pn = 
Pt = 


1 
1 
1 


19898 
20285.1 
20285.1 


24981.5 












Harigaya 


Pn = 


1 


14662 


20669.7 


25840 








(9,0) 


II 


Pt = 


1 


11263.1 


13193.4 


18246.8 


27226.6 








Toy 


Pn = 


1 


10989.6 


13350 


13974.8 


19703.3 


24117.3 


25765.3 




M 


Pt = 


1 


11123.5 


13318.4 


18508.1 


22206.5 


23535.6 


25611.1 



Table 1. Calculated optical transition frequencies in both models 
for open ended SWCNTs and SWCNTs with periodic conditions. 
For the Pn = 1 hnes /3s: = 1 and /3t = and for the hnes /3t = 1 
we have Ps = Pn = 1- 



Both models give the same band gap for (4,0) nanotube. It agrees well with the 
gap value given in the literature. For (9,0) the results of both models are the same. 
When the toy model is used one gets many more vHSs. For (9,0) we obtained an 
extra singularity in the vicinity of Fermi energy. The calculated next three singu- 
larities agree with those given in the literature. Wc could not find any data in the 
literature for the correlation of the other singularities. 
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